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Abstract. We define a new algebra, which can formally be considered as a CV deformed 
su(2) Lie algebra. Then, we present a one-dimensional quantum oscillator model, of which 
the wavefunctions of even and odd states are expressed by Krawtchouk polynomials with fixed 
p = 1/2, K 2 n{k', 1/2, 2j) and K^nik — 1; 1/2, 2/ — 2). The dynamical symmetry of the model 
is the newly introduced su(2)cv algebra. The model itself gives rise to a finite and discrete 
spectrum for all physical operators (such as position and momentum). Among the set of finite 
oscillator models it is unique in the sense that any specific limit reducing it to a known oscillator 
models does not exist. 


1. Introduction 

Krawtchouk polynomials are the simplest finite-discrete polynomials of the Askey scheme of 
orthogonal polynomials. They play an essential role in the study of problems coming from 
different branches of modern science. Image analysis based on their moments [T], exactly solvable 
birth and death processes [2], perfect qubit transfer in a spin chain with nearest-neighbour 
interaction [3] are only a few of such problems solved by employing Krawtchouk polynomials. 
Also in the “discretisation” of the quantum-mechanical formulation of the harmonic oscillator 
problem they play an important role. The non-relativistic explicit formulation of this problem 
in the infinite continuous configuration space is well known within the canonical approach. The 
wavefunctions of stationary states are expressed in terms of the Hermite polynomials both in 
momentum and position representation [3]. Under the assumption that the quantum world 
is quite different from the classical one, various approaches and methods allow to construct a 
number of exactly solvable mathematical models of the quantum harmonic oscillator, where 
the wavefunctions are expressed in terms of other known orthogonal polynomials. Also the 
Krawtchouk polynomials allow one to construct a number of interesting quantum harmonic 
oscillator models in finite-discrete configuration space. At present, the known models using 
Krawtchouk polynomials are the su(2) model [SHIj and the 5l(2|l) model [S] for the finite 
quantum harmonic oscillator. In the case of 5u(2) dynamical symmetry, the wavefunctions of the 
model are defined through the special case of Krawtchouk polynomials with a fixed parameter 
p = 1/2. For the ‘supersymmetric’ model with s[(2|l) dynamical symmetry, the wavefunctions 


are expressed by means of general Krawtchouk polynomials with parameter 0 < p < 1. In 
the present work, we present a new oscillator model with wavefunctions of even and odd states 
also expressed by Krawtchouk polynomials with fixed p = 1/2, but in a very different way: 
the polynomials that appear are K 2 nik', 1/2, 2j) and K 2 n{k — 1; 1/2, 2j — 2) respectively. The 
dynamical symmetry is formally called a CV deformation of the 5u(2) Lie algebra. 

The current paper is structured as follows; in section 2, the CV deformation of the su(2) Lie 
algebra and its representations are introduced. Then, this algebra is used for the construction of 
a model of the finite-discrete quantum harmonic oscillator. In section 4, we discuss explicit 
expression of the wavefunctions of the constructed model both in momentum and position 
representations. Finally, section 5 covers a brief summary and discussion of the obtained results. 


2. The algebra su(2) under the CV deformed symmetry 

Our starting point is the definition of the formal algebra su(2) under the CV deformed symmetry 
by its basis elements Jq, J+ and J_. 

Definition 1 The algebra 5u{2)c-p is a unital algebra with basis elements Jq, J_, C and V 
subject to the following relations: 

• The operator C commutes with all basis elements. 

• V is a parity operator satisfying V^ = 1 and 


[V,Jo]=VJo-JoV = 0, {V,J±}=VJ± + J±V = 0. (1) 

• The su(2) commutation relations are CV deformed as follows: 

[Jo, J±] = :^J±, (2) 

[J+,J_] = 2Jo(CiP-l). (3) 


The next step is to describe representations of this algebra, by defining the actions of the 
basis elements Jq, J+ and J_ in such a way that they are compatible with the commutation 
relations ([2]) and ([3]) . It is easy to verify that an action exists on the usual standard basis vectors 
\j, m) (with j being an integer and m = —j, —j -|- 1 ,... , +j), which is reminiscent of but different 
from the common su(2) algebra action. This is given by: 


Jo\j, m) 
J+\j, m) 


J-\j, m) 


= m \j,m), 


V(j - ^)(j - m-l) \ j,m + 1), 
V (i + m){j + m + 1) \ j,m + 1), 


\/(j + m){j + m- 1) \ j,m- 1), 
VU - '^)ij -m + l) \ j,m- 1), 


if j + m is even; 
if J + m is odd, 

if J + m is even; 
if j + m is odd. 


( 4 ) 

( 5 ) 

( 6 ) 


Note that j must be integer (and not half-integer) in order to be a representation (one way to 
see it is that one should still have J+\j,j) = 0 and J-\j, — j) = 0). The algebra representation 
is completed by the action of the extra operators C and V: 


C\j,m) = 2j\j,m), (7) 

V\j,m) = {-iy+^\j,m)- ( 8 ) 

Despite the fact that this algebra is at first sight very close to su(2), its behaviour is quite 
different. For example, the representations matrices of Jq, J± in the su(2)c-p algebra never 
coincide with the representation matrices of these elements in the su(2) algebra, except for the 
trivial representation j = 0. Note also that there is no “deformation parameter” that takes the 
su(2)c'P algebra to the su(2) algebra for some limit. 






3. A one-dimensional oscillator model based on the algebra su(2)cp 
We now consider an oscillator model based on the algebra introduced, and choose the 
Hamiltonian, position and momentum operators in a way that is similar to the non-deformed 
su(2) one-dimensional oscillator case [6l[7], i.e. 

H = JQ+j + -., q=-[J+ + J-), p=-(J+ —J_). (9) 

One can easily check that these three operators satisfy the Heisenberg equations: 


HA 


= -ip, 


HA 


= 19- 


Then, one can also observe that 


H\j,m) = {m + j + -)\j,m). 


So the spectrum of H is equidistant but finite: 

1 


n + - (n = 0, l,...,2j). 


( 10 ) 


( 11 ) 


( 12 ) 


From the actions ([S]) and ([^ one obtains for the position operator q that 


2q \j, m) = y(7+^rjy(JTmr^ |j, m - 1) + I j, m + l) , (13) 

when j + m = 2n is even, and 

2q \j, m) = |j, m - 1) + V U + m) {j + m + 1) |j, m + 1) , (14) 


when j -|- m = 2n -|- 1 is odd. One can do the same computations for the operator 2ip. In matrix 
form, both position and momentum operators 2q and 2ip take the following form: 


and 


2q = 


2ip = 


with matrix elements 


/ 0 Mq 0 

Mo 0 Ml 

0 Ml 0 


\ 0 0 


M2j-l 


f 0 Mo 0 

-Mo 0 Ml 


0 -Ml 0 


\ 0 


0 \ 
0 


M2j-1 

0 / 


-M2j-l 


= M'', 


0 \ 
0 


M2j-l 

0 / 


= MP, 


Mk = 


j ^/k {k + 1), if k is odd; 

I y/{2j — k) (2j — k — 1), if k is even. 


(15) 


(16) 


(17) 


The purpose is now to study the eigenvalues of the position (and momentum) operator in this 
model. In other words: we need to find eigenvalues and eigenvectors of these matrices and 
















M^. It turns out that these eigenvectors can be constructed in terms of symmetric Krawtchouk 
polynomials. In general, Krawtchouk polynomials Kn {x‘,p, N) of degree n {n = 0,1,..., iV) 
in the variable x, with parameter 0 < p < 1, are defined in terms of the 2 F 1 hypergeometric 
function of argument ^ as follows [9]: 


Kn{x-p,N)= 

Their orthogonality holds for discrete values of x, and is given by: 

N 

{x,N) Knix-p,N) Kn' {x;p,N) = h{n,N)6ny, (19) 


where, 


w (x, N) 


h{n,N) 


(x = 0,l,...,iV), 

(-l)"n! A-pV 
i-N)n [ P ) ' 


( 20 ) 

( 21 ) 


Here, {a)k = a{a + 1) • • • (a + A: — 1) are Pochhammer symbols. We can introduce the following 
orthonormal Krawtchouk functions: 


kn {x]p,N) = Y {x;p,N) . (22) 

Let us also mention that for p = 1/2, the corresponding Krawtchouk polynomials are called 
“symmetric”, since 

K„(x;i,iV) = (-l)-K„(iV-x;i,A^). (23) 

In previous examples, the construction of a quantum oscillator model in term of orthogonal 
polynomials is based on the differential (or difference) equation, the solution of which is 
expressed in terms of that polynomial. In this sense, the 5u(2) one-dimensional Krawtchouk 
oscillator model [SIE] is based on the difference equation for Krawtchouk polynomials m 
( 9 . 11 . 5 )], whereas the s[(2|l) one-dimensional Krawtchouk oscillator model [5] is based on 
the forward and backward shift operator equations for them (9.11.6), (9.11.8)]. For our 
current construction, the essential point here is a set of new difference equations for symmetric 
Krawtchouk polynomials {p = ^). These new relations involve a pair of symmetric Krawtchouk 
polynomials shifted in variables 2x or 2(x -|- 1) and parameters 2(j — 1) and 2j. 

Proposition 2 The symmetric Krawtchouk polynomials satisfy the following difference 
equations: 

j (2j - 1) Kj+n (2 (x -b 1); 2j) = - (x -b 1) (2x -b 1) Kj+n-i (2x; 2 (j - 1)) 

+ (j - X - 1) (2j - 2x - 3) Kj+n-i (2 (x + 1); i, 2 (j - 1)) ; (24) 

2 (j + n) (j - n) Kj+n-i (2x; i, 2 {j - 1)) = j (2j - 1) Kj+n (2x; i, 2j) 

- j {2j - 1) Kj+n (2 (x + 1); i, 2j) . (25) 

In these equations, the integer n is arbitrary but should he taken appropriately (e.g. for the first 
Krawtchouk polynomial, the degree j + n should belong to {0,1,..., 2j}). 






Proof. The validity of both equations is related to a higher level hypergeometric series than 
2 F 1 , namely 3 F 2 . For this purpose, we can use the following relation Uni (39)]: 


2 F 1 


f — 2 x, —j — n 

V -2i 




3 F 2 


/-n,n,-x \ 

V 1/2,-j ’ y ■ 


(26) 


Using this equation for the 2 F 1 expressions in (I24|) leads to the following identity to verify: 


j 3 F 2 


/ -n,n,-x - 1 \ 

V 1/2,-i ’ J 

-{x- j + 1 ) 3 F 2 ( 


= (a: + 1) 3 F 2 


/ -n,n,-x \ 

Vl/2,-J + l’ J 


—n, n, —X — 1 \ 

1/2,-j + 1 ’ )■ 


(27) 


Now, to prove the correctness of (1271) it is sufficient to simplify the Pochhammer symbols in 
the expansion of the 3 F 2 ’s in the right hand side, collect terms, and the left hand side appears 
automatically. 

The proof of (1251) is similar. In this case, using (|26p leads to the following identity: 


j'^ — in? 


3 F 2 


-n,n,-x A 1 f-n,n,-x-l 


1 . 


-{x-j + ^) 3F2 


—n, n, —X 

1 / 2 ,-j ’ 


(28) 


The verification of this again depends on elementary manipulations of Pochhammer symbols. □ 
Now, our aim is to construct eigenvalues and eigenvectors of M”? (1151) . This is rather technical, 
but is essentially based on (IMl) and (l25l) . For this purpose, we construct the matrix U in the 
following proposition. 

Proposition 3 Let AT? = 2q be the tridiagonal {2j + 1) x (2j + 1) -matrix (fT5]) and let 
U = {Uki)o<k,i< 2 j the (2j + 1) X {2j + l)-matrix with matrix elements as follows. For an 
even row index 2i (i = 0,1,..., j), let 

U2i,j-k = U2i,j+k = {-ly2^-^K2^{k■,\,2j), (29) 

U2i,o = U2i,2j = {-iy2-^/^K2i{j; i2i); 

C/2i,, = (-l)*2-l+i/2i^2,(oA,2j). 


For an odd row index 2i + 1 (i = 0,1,..., j — 1), let 

U2^+lJ-k = -U2i+i,j+k = {-iy+^2^-^ K2i{k - 1; 1, 2j -2), kG j-1} (30) 

U2i+i,o = -U2i+i,2j = (-l)'+'2-i/2^2i(j - 1; i2j - 2); 

U2i+l,j = 0 - 


Then U is an orthogonal matrix: 

UU^ = U^U = I. (31) 

Furthermore, the columns of U are the eigenveetors of , i.e. 


= UOy (32) 

where D'^ = diag(eo, ei,..., e 2 j) is a diagonal matrix eontaining the eigenvalues of TP?.- 


Efc = -2^/{j - k){j + k), ej = 0, e 2 j-k = 2^/{j - k){j + k), {k = 0,..., j - 1). (33) 






Proof. The proof of this proposition is quite similar to [101 Proposition 3]. The matrix 
relations (j32p follow directly from the two difference equations (12411 and ([25]). The orthogonality 
of the matrix U is slightly tricky, and does not follow (as one would expect) directly from 
the orthogonality of symmetric Krawtchouk polynomials. Instead, just as in the proof of 
Proposition 2 we need to go to higher level hyper geometric series. Reexpressing the 2^1 s in 
terms of 3 T 2 ’s by means of (12611 . the corresponding 3 T 2 ’s can be related to Hahn polynomials [9] 
with parameters a = jS = —1/2. The orthogonality of these Hahn polynomials leads to: 



2N)Kn'{2x-\aN) 


^n,n' 1 if 71 ^ A^, 

25n,n'-, ifra = A^, 


(34) 


for n, n' € {0,1,... , N}. The extra factor 2 for n = n' = N appears because for that case one 
has to take a limit of the squared norm of the Hahn polynomial, so one gets another factor than 
in the other cases where one can just substitute a = (3 = —1/2. Equation (I34|l is the underlying 
equation proving the orthogonality of the columns of U. The extra factor 2 is also responsible 
for the extra factors 2 ^/^ when k = j. □ 

This proposition gives us essentially the eigenvalues q of the position operator q. There are 
2j + 1 distinct eigenvalues and they are given by: 


- i, -Vf -1,... - - (j - l)^ 0, - {j - if, 1, j, ( 35 ) 

or in general: 

Q±{j-k) = fc = 0, l,2,...,j. (36) 

In Figure 1, we present a plot of the spectrum of the position operator for the su(2), s[(2|l) 
and su(2)c-p oscillator models. One observes that the spectrum of the position operator is quite 
different compared to the finite-discrete oscillator models already known. The main difference 
is the distribution of positions. For the su(2) oscillator model, the positions are equidistant 
{Q±k^'^ = 3zk). The distribution of positions for the s[(2|l) model is not equidistant, due to 

~ ±\/fc. In that case, the interval covering the positions is smaller than the corresponding 
one for the su(2) case. Comparing the 5u{2)c'p case with two others, we observe that some of 
its properties overlap with su(2) and others with sl(2|l). The interval covering the positions is 
the same for the 5u(2) and the su(2)cp model: ■ Both the s[(2|l) and su(2)c-p 

oscillator models have a non-equidistant distribution of the position eigenvalues. The main 
feature of the 5u{2)c'p oscillator model is that the position values are concentrated near the 
border of the inteval [—j, +j], whereas near the middle of the interval the distribution is sparse. 


4. Position and momentum wavefunctions 

The position wavefunctions of the su(2)cp finite-discrete oscillator model are defined in a 
similar way as in [101 • Their explicit expression is computed through the overlap between 
the ^-eigenvectors \j,qj-k) and the ^/-eigenvectors \j,m). Let us denote them by 
where m = —j, —j -|- 1,..., +j, and where q assumes one of the discrete values of qk with 
{k = —j, —j + l,...,+j)- Simple computations show that 

'kj+m((?j—fc) {j:'^\j:Qj—k) f/j+mj+fc) (3^) 

where Uj+m,j+k are the matrix elements of M'^. Then, one can find that for the even case 
j + m = 2n and positive values of position, the wavefunctions have the following expression: 


^2n(%) = (-1)''2^ ^R2n(fc; ^,2j), 77 = 0, 1, . . . , j, A: = l,...,j-1 


(38) 
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Figure 1. Spectrum of the position operator for j = 5, (a) in the case of the su(2) model, (b) in 
the case of the s[(2|l) model and (c) in the case of the su(2)c'P model. 


whereas, for the odd case j + m = 2n + 1 and positive values of position, we have 

^ 2 n+i{qk) = {-lT 2 ^-^k 2 n{k - 1; i, 2j - 2). (39) 

One can extend these computations and obtain similar expressions for zero and negative values 
of the positions. 

In Figure 2 we plot the discrete wavefunctions (|38p and (j39p for some large j-value, and for 
some n-values. Due to the unique property of the position eigenvalues to be located near the 
border of the interval [—+i], the behaviour of the wavefunctions also reflects this peculair 
distrubution. 

We can define momentum wavefunctions in a similar way. Their explicit expression 

is computed through the overlap between the p-eigenvectors \j-,Pj-k) and the iJ-eigenvectors 


|j, m). In this case, the simple relation between and shows that 

^j+miPj—k) (j) ^|j) Pj—fc) ^j+m,j+ki Pj—k — Qj—k^ ^k,l — 1 (^0) 

Then, in accordance with equations (j38p and (13911 . one finds that for the even case j + m = 2n 
and for positive values of momentum, the wavefunctions have the following expression: 

^2niPk) = ^2^~^k2n{k;l,2j), n = 0, l,...,j. A: = 1,..., j - 1, (41) 

whereas for the odd case j + m = 2re + 1 and positive values of momentum, we have 

^2n+l{Pk) = -2^-^k2n{k - V, ^ 2j - 2). (42) 

These computations can easily be extended for zero and negative values of the momentum. Both 
wavefunctions 'l'(gr) and 4>(p) are related by the so called CV symmetric Krawtchouk transform; 

j 

^{pi)= (43) 

k=-j 

where, /C = U'^V. The matrix elements JCk^i can be computed by using (IM|1 . 





Figure 2. Plots of the discrete wavefunctions ^n{Q) in the representation with j = 30. We plot 
the wavefunctions for n = 0, 1 , 2 , 60. 


5. Discussion 

Despite the fact that the non-relativistic one-dimensional quantum oscillator in the canonical 
approach has a nice explicit solution in terms of Hermite polynomials, there are still attempts 
to construct new oscillator models inspired by different starting points. One of these starting 
points is the assumption of discrete space, in which case one is looking for oscillator models with 
a discrete position spectrum (which could be hnite of infinite discrete). The most popular finite 
discrete oscillator model mm is based on the su(2) algebra, and uses Krawtchouk polynomials. 

In current paper, we have constructed a new model for the one-dimensional oscillator, with 
wavefunctions expressed in terms of the Krawtchouk polynomials with parameter p = 1/2 (the 
symmetric Krawtchouk polynomials). In order to construct a solvable oscillator model, we have 
first defined a new algebra, formally called the su(2) algebra under the CV deformation. The 
representations of this algebra were given and are similar to representations of 5u(2), but the 
action of the raising and lowering operators J+ and J_ differs considerably. As a consequence, 
when considering the oscillator model related to this algebra, the matrix representations of 
position and momentum operators is also different. In the standard \j, m) bases, the position 
operator is tridiagonal, and we have shown that its eigenvalues and eigenvectors involve 
Krawtchouk polynomials. The matrix relation for the eigenvalue problem is equivalent to certain 
finite difference equations for Krawtchouk polynomials. 

The resulting model is unique in the sense that it does not generalize any other existing 
oscillator model. Peculiar properties are the non-equidistant distribution of the position values 
as well as the behaviour of the wavefunctions around zero and near the border of the maximum 
and minimum eigenvalues. 
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